Analysis of Shear-Thickening in Physical Gel. 
A Stochastic Theory for Polymer Networks 
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A formula of steady shear viscosity is derived by introducing a model to describe the dynamics 
of physically cross-linked network (physical gel), and successfully analyzes the shear-thickening 
behavior observed in HEUR aqueous solutions by Jenkins, Sileibi and El-Aasser. We take into 
account the effects of looped chains at junctions which detach their one end from the junction as 
the shear rate increases due to the collisions with other chains. This process produces the weak and 
wide enhancement of the number of chains whose both ends stick to the separate junctions (active 
chains) leading to the weak increase in the shear viscosity. It is also shown that the nonlinear force 
sustained by active chains induces the strong and sharp enhancement of the steady shear viscosity. 
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I. INTRODUCTION 

One observes generally that the viscosity of poly- 
mer solution monotonically decreases when one increases 
shear rate. This effect, referred to as shear-thinning, re- 
sults from the disentanglement of polymer chains under 
the shear flow. However, some sort of polymer solutions 
exhibit an unusual behavior, called shear-thickening, i.e., 
the viscosity of the solution grows as the shear rate in- 
creases, attains a maximum and decreases at higher shear 
rates. 

The shear-thickening phenomenon is observed for poly- 
mers having a few segments capable of association, such 
as hydKtt^jJiCp,aolymers containing a few hydrophobic 
groups ,BffloEoEl and ionomersBcl In solution, these asso- 
ciative polymers construct a physically cross-linked net- 
work due to temporary connections among their associa- 
tive segments (physical gel). A typical example of such a 
physical gel is HEUR (hydrophobkalhj-mDdified ethoxy- 
lated urethane) aqueous solutionffiMtffl HEUR is a 
hydrophilic poly(ethylene oxide) containing hydrophobic 
hydrocarbon groups at its both ends. In water, the hy- 
drocarbon groups at their ends associate with each other 
due to hydrophobic interactions among them to construct 
a physically cross-linked network at relatively low poly- 
mer concentrations. Another example is a solution of 
telechelic ionomers in apolar solvent Jju such as a solution 
of a, w-Mg carboxylato-polyisoprene in decahydronaph- 
thaleneE 

In this paper, we study linear chains having associa- 
tive functional groups at both their ends, and analyze 
the steady shear viscosity of HEUR aqueous solutions o. 
served by Jenkins et alB representing shear-thickening.l 
We will consider three types of chains in the model, i.e., 
active chains with its both ends connecting to separate 
junctions, dangling chains with only one of its end stick- 
ing to a junction and loops with both ends connecting 
to a junction (Figure |l|). Floating (free) chains in sol- 
vent are not taken into account in the present approach. 
We introduce loops since HEUR in water has been ob- 



served to form flowerlike micelles comprised primarily of 
loops, as schematically depicted in Figure [l], for broad 
polymer concentrations above the onset of association^ 
A dangling chain can become either an active chain or 
a loop by connecting, respectively, its one free end to a 
junction with which its other end does not stick or has al- 
ready stuck. In other words, an active chain can become 
a loop via dangling state and vice verse. We take into 
account the breaking of a loop by detaching its one end 
from a junction due to the collisions with dangling chains 
or other loops as the shear rate increases. It will be seen 
that this process produces the weak and wide enhance- 
ment of the number of active chains, leading to the weak 
increase in the shear viscosity. The idea of increase in 
the number of active chains around the shear-thickening 
region has been suggested from an observation that the 
plateau modulus there is larger than that in the New- 
tonian region.Q It will be also shown that the nonlinear 
force of active chains described by the inverse Langevin 
function induces the strong and sharp enhancement of 
the steady shear viscosity. The effect of the nonlinear 
force has been discussed by others in diffeEent approaches 
both theoreticallyEH and experimentallytj. The success 
in the analysis of the shear viscosity in HEUR aqueous 
solutions by the present theoretical formula reveals that 
both the increase in the number of active chains and the 
nonlinear force of active chains are indispensable to un- 
derstand the shear-thickening behavior. 

Physical gel is a polymer network in which junctions 
can break and recombine due to thermal fluctuation, i.e., 
the junctions have the finite lifetime. The transient net- 
work theory, formed first by Green and Tobolsky,E-H is a 
powerful tool to investigate the rheological property of 
such gels. They extended the classical model of rubber 
elasticity so as to allow for the junctions which break and 
reform with constant rate independent of shear rate. Due 
to the assumptions of the constant breaking rate of junc- 
tions and of the Gaussian chains, they could not succeed, 
to derive any non-linear rheological behavior. Lodgelij 
generalized the Green- Tobolsky (GT) model by intro- 
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FIG. 1: Schematic representation of a network constructed 
by HEUR in aqueous solution.!] 



ducing the distribution in the lifetime of junctions (the 
inverse of breaking rate of junctions) which does not de- 
pend on the shear rate. Like the GT model, obtained 
shear viscosity did not depend on the shear rate, i.e., 
there was neither shear-thinning nor shear-thickening. 
On the other hand, YamamotoEj evolved the GT model 
by introducing the breakage rate of junctions which de- 
pends on the end-to-end length of a section of chain lying 
between two neighboring junctions (partial chain). Since 
the flow changes the end-to-end length of partial chains, 
the breakage rate must depend on the rate of flow. By as- 
suming that the junction breaks whenever the end-to-end 
length of the partial chain exceeds a certain critical value 
under the flow and that the polymer chain is Gaussian, 
he showed that shear-thinning appears. This is because 
the number of partial chains decreases as the shear rate 
grows. 

With the intention to analyze the physical gel com- 
prised of linear chains having functional grou-ps-at their 
each end, like HEUR, Tanaka and EdwardstZllia devel- 
oped the transient network model formed by Yamamoto 
by introducing two types of chains in the network, i.e., 
the active chain and the dangling chain. Similar to the 
Yamamoto model, they adopted the breakage rate de- 
pending on the end-to-end length of active chains. As- 
suming that the breakage rate is an increasing function of 
the end-to-end length of an active chain without a cutoff, 
and that the active chain is Gaussian with linear force, 
they obtained the shear-thinning behavior, but could not 
the shear-thickening. 

Wang0 extended the Tanaka-Edwards (TE) model by 
introducing free chains with both ends being not con- 
nected to any junctions in addition to active and dan- 
gling chains. Like Tanaka and Edwards, he adopted the 
increasing breakage rate of an active chain without intro- 
ducing a cutoff to its end-to-end length. He assumed that 
shear flow promotes creation of more free chains which 
collide with junctions resulting in the growth of the tran- 
sition rate from a free chain to a dangling chain. Hence, 
the number of dangling chains grows as the shear rate in- 



creases, leading to the enhancement of the number of ac- 
tive chains which induces shear-thickening. However, the 
adoption of this transition rate appears to be somewhat 
doubtful, since it may not be affected by the frequency 
of collisions of a free chain with junctions but affected 
by the,»umber of junctions surrounding the ends of free 
chains S3 

All the transient network models mentioned above 
have one common assumption, i.e., the relaxation time 
of dangling (and free) chains is much smaller than that 
of the network. Under the assumption, a dangling chain 
immediately relaxes to, say, the Gaussian chain after its 
connected-end detaches from a junction. Vaccaro and 
MarruccinJ considered the network model composed of 
active and dangling chains, however, in order to avoid 
the above assumption, they introduced a 'Fokkcr-Planck 
equation' for the distribution of the end-to-end vector 
of daagling chains. They adopted the nonlinear Warner 
forceEil for active chains, and introduced the breakage 
rate of active chains which increases with their end-to- 
end length. They also introduced the recombination rate 
linearly increasing with respect to the end-to-end length 
of dangling chains. For the relaxation time r of chains 
having the order of the lifetime (at equilibrium state) 
/3g of active chains, shear-thickening is show n to ap- 
pear. On the other hand, as seen in subsection VB, the 



relaxation timej of chains in the experiment we are go- 
ing to analyzed is estimated as r ~ 10~ 5 ~ 10~ 4 s, while 
Pq 1 ~ 1CU 2 ~ 10 _1 s. For such short relaxation time 
(r <C Pq 1 ), shear-thickening does not seem to appear 
within their model. Therefore, we need to propose other 
mechanism to explain the shear-thickening phenomenon 
even in the case t^Pq 1 . 

The transient network theory may be also applicable 
to the topologically entangled polymer system when the 
entanglement points are regarded as junctions since the 
entanglements formed by chains can be disentangled due 
to the thermal agitation of chains along the tube-like re- 
gion comprised of suipounding chains, i.e., reptation pro- 
posed by de Gennes.E3 Actually, Yamamoto evolved his 
model with the intention to treat polymer solutions at 
rather high concentration. However, the topological re- 
striction of chains is only effectivelv-taken into account in 
this treatment. Doi and EdwardsEa developed a theory 
for the concentrated polymer solutions or melts on the 
basis of the reptation model. They succeeded, for exam- 
ple, to derive shear-thinning behavior for a topologically 
entangled system comprised of linear chains. 

In the present paper, we will focus our attention on the 
system where the above effect of the topological entan- 
glement is week and negligible, since in the experiment 
we will analyze, the polymer concentration is at most 
on the order of the overlap threshold of polymer chains 
(see subsection VA). Junctions are ascribed to the as- 
sociation among functional groups at each end of linear 
chains. 

In section H, we introduce the basic equation deal- 
ing with the physically cross-linked network composed 
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of unentangled linear chains with associative functional 
groups at both their ends. The derived equation, rep- 
resenting the transition among each type of the chains, 
will be used throughout in this paper for the analysis of 



the steady shear viscosity of physical gel. In section II] 



the dynamics of each type of chains and the transition 
rates among them under the steady shear flow will be 
discussed. We see that the transition from a loop to a 
dangling chain is promoted by the shear flow. It plays an 
important role in the shear-thickening phenomena as well 
as the nonlinear force of active chains. In section |y|, we 
solve the basic equation introduced in section |l| in order 
to give a formula for the steady shear viscosity. In section 
[v| we analyze the experimentally observed shear viscos- 
ity of HEUR aqueous solutions^ by the obtained formula. 
It will be shown that the formula explains experiments 
quite well for the broad shear rate. Discussions will be 
devoted to section VI. 



II. BASIC EQUATION 
A. Assumptions 

We treat in this paper the system made up of linear 
chains which are uniformly distributed in space, with 
uniform molecular weight, carrying associative groups at 
both their ends. No topological entanglements are con- 
sidered. Here we list again the three types of chains rele- 
vant to our model (Figure |l|): 1) Active chain: The both 
ends of the chain are connected to two junctions sepa- 
rated in space. 2) Dangling chain: Only one of the two 
ends is connected to a junction. 3) Loop: Both ends are 
connected to a single junction. 

Suppose the situation in which an incompressible 
macro-deformation is added to the uniform system. We 
assume that the distribution of chains remains uniform 
even after the application of the deformation. Hereafter, 
we will pay our attention to those phenomena taking 
place in an arbitrary unit volume in the system. The 
number of chains n per unite volume is constant in time 
thanks to the incompressibility of the deformation, and 
in space due to the homogeneity of the system. 

The added deformation produces a velocity-field in the 
system causing a conformation change of chains. We rep- 
resent the conformation of a chain by means of its end- 
to-end vector. 



B. End-to-End Vector Space 

Now, we introduce a space, named the end-to-end vec- 
tor space, in which we will construct a basic equation 
for the end-to-end vector of chains. In this space, all 
chains which have the same end-to-end vector are repre- 
sented by one point called a representative point for the 
chains. Since our system consists of those chains which 
have the same contour length, say I, all representative 



points locate within the sphere of radius I. Note that the 
representative point have nothing to do with a position 
of chain in real space. Even if a chain moves around in 
real space, its representative point does not move unless 
the chain changes its end-to-end vector. 

Thermodynamical quantities related to a chain, such 
as tension, free energy and so on, are often given by func- 
tions of its end-to-end vector r. The end-to-end vector 
space provides us with the support of these functions. 



C. Equation for the Number of Chains 

1. Eulerian Description 

Let us introduce, in unit volume in real space, the num- 
ber (f> l (r, t)dv of z-chains with end-to-end vector r^r + dr 
at time t. The superscript i distinguishes the type of 
chains, i.e., i = a represents active chains, i = d dan- 
gling chains and i — I loops. Note that the total number 
4>(r, t)dr of chains with end-to-end vector r ~ r + dr in 
unit volume is given by 



(r,i)dr = ^0 l (r,i)dr. 



(2.1) 



We derive the equation for 4> z (r, t) in the end-to-end 
vector space for z-chains under the velocity-field r l (r,t) 
caused by the deformation added to the system. 

Let us consider a volume V 1 (t) , in the end-to-end vec- 
tor space for z-chains, which varies its boundary in accor- 
dance with r z (r,t). Since the points for z-chains in this 
volume are created or annihilated, we have the equation 
for number-conservation of z-chains 

( / dv 4> l {v,t) ) = - ]T / dr W ot (v,tW(v,t) 

+ J2 [ dr Wy(r,t)^'(r,t), (2.2) 

where Wjj(r, t)dt is the transition probability during dt 
at time t from j-chains to z-chains having the same end- 
to-end vector r (see Figure ||). The conservation of the 
transition probability Wij (r, t)dt reads 



J2Wij{r,t)dt = l. 



(2.3) 



The first term in the right-hand side of (2J2) represents 
the number of z-chains in V l (t) which convert to other 
types of chains in unit time, whereas the second term 
stands for the number of z)-chains in V z (t) which 
convert to z-chains in unit time. That is, the first and 
the second terms represent annihilation and creation of 
z-chains, respectively. As V l (t) is arbitrary, we obtain 
the equation for (j) l (r,t) from (|2.2[) 



r l (r,t)^(r,i) 
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= - £ WflWfct) + E ^ WM)- (2.4) 

When the ioiaZ number of chains does not change, 
(r, t) satisfies an equation 



j^(r,t) + ^- (r(r,t)^(r,«)J =0. (2.5) 



The requirement that the sum mat ion of (2.4) with re- 
spect to i should coincide with (2.5) gives us the relation 
between r(r, t) and r l (r, t) in the form 



r(r,t)(f>(r,t) ^r'fr,^ (r,(). 




Y Wji(li,t) 

Wij(ri,t)A _ T 0J(r,t) 




FIG. 2: Schematic explanation of transitions between i- 
chains and j-chains having the same end-to-end vector ri 
at time t (Eulerian description). The upper plane (actually 
three dimensional) stands for the end-to-end vector space for 
i-chains and the lower for ji-chains. A solid line in the end-to- 
end vector space for a-chains (a — i,j) represents a streamline 
for a-chains passing through a point ri at t. The streamline 
gives the velocity r Q (r, t) as its tangent line at r. A represen- 
tative point for a-chains moves in accordance with r a (r, t) in 
each end-to-end vector space as described by (2.4). 



2. Lagrangian Description 

The basic equation (|]J) of the present paper is the one 
in the Eulerian description. Since, hereafter, we will an- 
alyze mainly in the Lagrangian description, we introduce 
the quantity 

4>\v l (t; 4, t ), t) =Jdr 5(r - v l (t; r , i o ))0 4 (r, t), (2.7) 

in the description where r 8 (i;r ,io) is a position of the 
representative point belonging to i-chains at t, which was 



located at r at the initial time to- rI (i; r cb^o) is a solu- 
tion of the equation of motion 



dr l (t) 
dt 



rV(t),t), 



(2.8) 



with the in itia l condition r*(to; r o; ^o) = r o- The equation 
of motion (2JS) describes the dynamics of the end-to-end 
vector in real space, which corresponds to the movement 
of a representative point under the influence of the flow 
r z (r, t) in the end-to-end vector space. 

Now, we derive the equation for </> l (r*(i; r ,t ),t). Dif- 
ferentiating (2.7) with respect to t, we have 



(2.6) d^(r*(t),t) 



dJ*(t,t ] 




J l (t,t ) dt 

(r,i) d 



dr 



where J l (t,to) is the Jacobian defined by 

drUt) 



^(r\t),t) 



J r (t,to) 



dr, 



00 



(2.10) 



For simplicity, we are omi tting r and to in r z (t;r ,to)- 
By putting th e eq uation (2.4) within the Eulerian de- 
scription into ( |2.9| ) and by multiplying J l (t,to) on both 
sides, (2.S) reduces to 



= -E W 3l {v\t),t)J\t,toW{v\t),t) 

+ £ WijiT^QjifatoWPfat). (2.11) 

A transition to the i-chain with end-to-end vector r l (t) 
(or to the region dr l {t) = J l (t,t )dr in the end-to-end 
vector for i-chains) can occur only from the j-chains 
U i) whose end-to-end vector r J (i) is equal to r l {t) 
(or from the end-to-end vector space for j -chains whose 
region dr>(t) — Ji (t,to)dr J is equal to dr'(t)). Therefore, 
the arguments of Wjj a nd </>> in the second term of the 
right-hand side of ( [2.11 ) are r l (t) instead of r? (t) (or the 
superscript of the Jacobian J l (t, to) is i instead of j) (see 
Figure ||) . 



Integrating (2.11) from time to to t, we have 
r(t,toW(r l (t;r l ,t ),t) 



— exp 



- I dt ' E W^v^t'^to)^') 

Jtn ,a 



^(4*0) 



/ di'cxp - / A" E ^(r l (i";r ,i ),t") 
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(2.4) becomes 



:drj(t)=JJ(t,to)dr J 



FIG. 3: Schematic explanation of the transition from j- 
chains to i-chains (Lagrangian description) . A solid line in the 
end-to-end vector space for i-chains (upper plane) expresses a 
trajectory of a representative point for i-chains which was at 
r at time to- On the other hand, a solid line in the end-to-end 
vector space for j-chains (lower plane) represents a trajectory 
of a representative point for j-chains passing through a point 
r 3 (t) which is equal to r l (t) at the moment t. The dark-gray 
region at r l (t) in the upper plane represents a domain dr 1 (t) = 
J 1 (t , to ) dr l Q whereas the light-gray region around r J (t) ( = r 8 (t) ) 
in the lower plain is dr 3 (t) = J J (t,to)dr J which is equal to 
dr l (t). 



x Y, ^ ii (r i (i / ;r^to),t / )J i (i / ,io)^'(r i (t';4,to),0- 

(2.12) 



The first term in the right-hand side of ( 2.12 ) expresses 
the number of j-chains with end-to-end vector r l (t) ~ 
r l (t) + dr 1 (t) which were i-chain having the end-to-end 
vector r <~ r +<irQ at time to, and remain i-chain un- 
til time t. This term vanishes at t — *■ oo in most cases. 
The second term represents the number of i-chains with 
r l (t) ~ r l (t) + dr 1 (t) which became i-chain at t' (> to), 
and remain i-chain until time t. 



= -J2 ^(r,t)4(r,t) + Y Wy(r,i)^(r,f), 

(2.13) 

and by replacing (p l (r l (t),t) and r J (t) with </>^(r^(t),i) 
and r^(i) respectively, (2.12) becomes 



r(t,toW f (r)(t;rlt ),t) 



exp 



I dt' Y, W^(r)(t';r ,t ),O 



rt 



+ f dt'exp - f dt" Y W ji (r)(t";rlt ),t") 

J to [ Jf j{lti) 

x Y W«(r/(*';'o J *o),0^(* / .*o)^(r>(t';r*,to) J 0» 

(2.14) 

where r l Jt; r , to) is a solution of the Langevin equation 



dr^t) 
dt 



r i (r i (t),t;K(t))=v}(r i (t),t), (2.15) 



with the initial condition ry(to; r , to) = r . An equation 
( 2.13 ) corresponds to the 'stochastic Liouville equation' 
in the end-to-end vector space for i-chains.t3 In the fol- 
lowing, we will denote by ((• ■ ■)) the average with respect 
to the stochastic process. 



E. Remarks 

In the present model, features of each type of chains 
are characterized through rl(r,t) and Wij(r, t). We will 
give functional forms of ij(r, t) and Wy(r, t) in the next 
section when the system is in the incompressible velocity 
gradient k(t), especially, in the steady shear flow. 



D. Stochastic Process 

When the dynamics of i-chains is given by a 
stochastic process, say R(t), r'(r, t) becomes stochastic 
r 8 (r, £; R(t)). In this case, 0''(r t) becomes also stochas- 
tic <fi l (r, t; R(f)) through (2.4). Hereafter, we denote 
r*(r,t;R(t)) as r}(r,t) and <f> l (r, t; R(t)) as </>}(r,t) for 
simplicity. The equations for </> l (r, t) in the previous sub- 
section hold for 0j (r, t) as well, if the products between 
two stochastic quantities, such as r^(r, t)d&(r, t), are re- 
garded as the Stratonovich-type products. E-3 By replacing 
4> l (r,t) and r l (r,t) with <^(r, t) and r\(r,t) respectively, 



III. CHARACTERIZATION 

A. Tension in a Chain 

When the distance between two ends of a flexible chain 
is finite, say r, there appears an effective attractive force 
f(r) between these ends, due to the thermal motion of 
segments within the chain. We adopt, for the attractive 
force, the one obtained within the random-flight model: 



f(r) 



a ^ r 



(3.1) 
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where L^ 1 is the inverse function of the Langevin func- 
tion defined by L(x) = cothx — l/x. Here and hereafter, 
TV is the number of statistical segments in a chain, a the 
length of a segment, T temperature of the solution in 
which chains are immersed and ks the Boltzmann con- 
stant. Note that in the random- flight model, the finitc- 
ness of the contour length I = Na of a chain is taken into 
account seriously. 



B. Dynamics in End-to-End Vector Space 

1. Active Chain 

We assume that all junctions move in accordance with 
the velocity gradient k(t), and that the effects of the fluc- 
tuation of junctions are week and negligible (affine de- 
formation assumption). Since the end points x"(t) and 
x 2 (*) °f an active chain in real space are sticking to differ- 
ent junctions by definition, the end-to-end vector r a (t), 
defined by r a (t) =x"(t) — xS^t), has finite value whose 
time evolution is given by 



dr a (t) 
dt 



= k{t)v a (t). 



(3.2) 



The equation (3.2) is equivalent to the velocity-field 



r a {r,t) = k(t)r 



(3.3) 



in the end-to-end vector space for active chains. It means 
that if a representative point at r belongs to active chains 
at t, it flows in accordance with r a (r, t). Integrating ( p 7 ^ ) 
from time to to t, we have 



r a (t;rg,i ) = Texp 



dt'k{t') 



rg = A(*,to)rg,(3.4) 



where T is the time-ordered operator. In the case of 
steady flows, since k is independent of time, X(t; t') be- 
comes a function of the time interval t — t' , i.e., A(i; t') = 
exp[(t-t')k] = X(t-t'). 



2. Dangling Chain 

One end of a dangling chain connects to a junction (we 
call it connected end), whereas the other end does not 
(we call it free end) . We treat the free end as a Brownian 
particle, i.e., the free end changes its position randomly 
by the random force R(t) acting on it. Let xf stands for 
the position of the free end of a dangling chain, and x. d 
the position of its connected end in real space. Dynamics 
of xf is determined by the Langevin equation 



,N 



dt 



k(t)4(t) = f (x?(t) - 4(t)) + R(t), (3.5) 



where we neglected the inertial term (oc d 2 xf(t)/dt 2 ) 
since the relaxation time for the velocity of typical poly- 
mer chains are usually much smaller than the time region 
of our interest Ej (n is the friction constant between a 
chain and a solvent, which is assumed to be Cn — -WCi 
where £i is the friction constant per segment. £i is given 
by Stokes's law: d = 6irr] s a with rj s being the viscosity of 
the solvent. The random force R(t) in (3.5) is assumed 
to be the Gaussian white process, i.e., its average and 
variance are given by 

((R a (t))) = 0, (3.6) 
((R a (t)Rp(s))) = 2( N k B T5 a!3 5(t - s), (3.7) 

respectively (a, (3 = x,y, z). The dynamics of Xj is deter- 
mined in accordance with k(t): 



dt 



k(t)?4(t). 



(3.8) 



Making use of ( |3.5| ) and ( |3.8| ), we obtain the equation 
for the time evolution of the end-to-end vector r d (t) = 
xf(t) — x^(i) in the form 



dr d (t) 
dt 



= n(t)r d (t) + l-^(r d (t))+K(t) 



(3.9) 



The equation (|3.9|) is equivalent to the velocity-field 



1 

0^ 



rj(r,t)=k(t)r+^[f(r)+R{t) ) (3.10) 



in the end-to-end vector space for dangling chains. It 
means that if a representative point at r belongs to dan- 
gling chains at t, it flows in accordance with r*(r, t). 

When the Hookean approximation is adopted for the 
tension f(r), i.e., 



f(r)=-^r (=f&(r)), K 



N 



the equation of motion (3.9) and the corresponding 
velocity-field (3.1C) become 



and 



dr d (t) 
dt 



rj(r,t) 



( £(t) - ^ W) + -^R(t) (3.12) 

y (jv j (,n 



k(t)-^-)v+ J-R(t), 



(3.13) 



respectively. The equation of motion ( 3.12 ) is analyti- 
cally integrated from to to t as 

rj^4,t )=e-^- t ^X(t,to)r d 

^' e - (t -*' )/2T A(i,t')R(i'), (3-14) 



CiVJt 
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where r is the relaxation time of a dangling chain defined 
by 



Cn CiN 2 a 2 



2K 



N 



6k B T 



(3.15) 



The Hookean approximation (3.11) is valid when the end- 
to-end length of a dangling chain is short enough com- 
pared with its contour le ngth. This situation is realized 
when the relaxation time (3.15) for the end-to-end vector 
of a dangling chain is much shorter than the characteris- 
tic time of the deformation. 



3. Loop 

The end-to-end distance of a loop is always zero, i.e., 

r l {t;r l Q = Q,to) = 0, (for ajiy t(> t )). (3.16) 

The equation ( |3.16| ) is equivalent to the velocity-field 

r'(r = 0,t) = 0, (for any t). (3.17) 

in the end-to-end vector space for loops. It means that 
if a representative point at r = belongs to loops at 
t, it does not flow. Note that no inform ation abo ut th e 
conformation of a chain is obtained from (3.16) or ( \i . 1 7 ) . 

Now, we estimate the dimension of a loop in a steady 
shear flow. Let x(n, t) be a position of the nth seg- 
ment (n = 1 , ■ • • , N) of a chain in real space at time 
t. The Langevin equation for segments is written, for 
?i = 2, ••• , N—l, as 



t / dx(n, t) „ . . . , 



ff (x(n, 

f)-x(n + l, + (*(«,*) 

— x(n 

+ R(M) 
<9 2 x(M) 



dn 2 



R(n,t) 



M)) 
(3.18) 

(3.19) 



where, in ( |3.1g ), we substituted (3.11) for the force fjy =1 
between neighboring two segments along the chain, and 
took the continuous limit for n. The random force 
R(n, t), acting on the segment at x(n, t), is assumed to be 
the Gaussian white process, i.e., its average and variance 
are given by 

{(R a (n,t)))=0, (3.20) 
((R a (n, t)Rp{m, s)))=2Q 1 k B T8 a f ) 8{n - m)S(t - a), (3.21) 

respectively. Since one end of a loop connects with its 
another end, the boundary condition at n — and iV = 
should be 



x(n = 0, t) = x(n= N, t) = 0, 



(3.22) 



where we fixed the connected point at the origin without 
loss of generality. The property of a loop is characterized 
through ( fT2^ ). 



Making use of (3.19) with (3.22), we can get the mean- 
square radius of gyration s a sp of a loop defined by 



S a Sf3 



N 



— / dn(([x a (n,t) - XGa(n,t)][xp(n,i) - xgp(ti, 
Jo 



(3.23) 



under a deformation. Here, X(j(t) = N^ 1 ^ dnx(n,t) is 
the center of mass of a loop. When a steady shear defor- 
mation with shear rate 7, described by the deformation 
tensor 



1 j(t-t') 
\{t,t') = I 1 
1 

or the velocity gradient tensor 

7 0' 
h =|000 




(3.24) 



(3.25) 



is added, (3.23) become (see appendix^) 



Na 2 



1 



13n 4 



2 _ 2 

S y — S z 



36 \ 10080 
Na 2 



{in? 



Na 2 
~12 



0. 



T\ is the relaxation time of a loop given by 

Ci^V 2 a 2 



2>Tt 2 k B T' 



(3.26) 

(3.27) 
(3.28) 

(3.29) 



which is essentially equivalent to r defined by (3.15) 



Transition Probabilities 



The j-chain becomes the i-chain with the probability 
Wij during unit time. We suppose that the transition 
rates are affected by how many times the j-chain 
collides with other chains within its lifetime under the 
influence of a flow. Since, in most cases, the frequency 
of the collision depends on the velocity gradient k(t) of 
a flow, Wij is a function of t through k(t) as well as the 
end-to-end vector r of the j-chain. However, at steady 
states in steady flows, Wij becomes a function of a quan- 
tity characterizing the flow, e.g., the shear rate 7 for the 
steady sh ear fl ow. In this subsection, for the shear flow 
given by ( 3.25| ), we introduce the dependence of Wij on 
7 and r. 

In the following, for simplicity, we denote the transition 
Wda from an active to a dangling chain by (3, W a d from 
a dangling to an active chain by p, Wm from a dangling 
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chain to a loop by v and Wdi from a loop to a dangling 
chain by u. We can put W\ a — W a i —0 between an active 
chain and a loop, since the primitive transitions should 
occur between the dangling and other states as stated in 
section QL 



1. Transition between Active and Dangling Chains 

a. From Active to Dangling Chains An active chain 
becomes a dangling chain with the transition rate (3 by 
detaching its one end from a junction. In our model, (3 
depends scarcely on 7. The reason is as follows. All active 
chains alter their end-to-end vector affinely, therefore, 
they hardly collide with each other. As for the dangling 
chain (or the loop), since it sticks to a certain junction, 
its center of mass moves in accordance with k. That is, 
on a certain xz-plane (i.e., the plane on which the shear 
flow has the same velocity) , the relative velocity between 
the active chain and the dangling chain (or the loop) is 
very slow. This is the reason why the frequency of the 
collision between the active chain and the dangling chain 
(or the loop) is small. 

We assume that the one end of an active chain dis- 
sociates from the network definitely when its end-to-end 
distance becomes longer than r*(< Na). We further as- 
sume that the breakage rate /3(r) for r < r* is constant 
and has value /3o- The breakage rate Sir) under these 
assumptions can be written in the formLS 



/3(r) = 



00 (r<r*) 
00 (r > r*) 



(3.30) 



b. From Dangling to Active Chains A dangling 
chain becomes an active chain with the transition rate 
p by attaching its free end to a certain junction with 
which its other end of the chain does not connect. Al- 
though the dangling chain collides with other dangling 
chains and loops, there seems to be no reason for the col- 
lision to affect the connecting process of the free end to 
the junction. Therefore, we assume that p is independent 
of 7. As for the r dependence, we can regard p as a con- 
stant effectively, since the connecting process does not 
seem to be affected strongly by the end-to-end vector. 



2. Transition between Dangling Chains and Loops 

a. From Dangling Chains to Loops A dangling chain 
becomes a loop with the transition rate v by attaching 
its free end to the junction with which its other end has 
already stuck. Since the connecting process is not influ- 
enced by the collision with other chains as stated above, 
v(r) can be assumed to be independent of 7. 

Since the process takes place only when its end-to-end 
length is 0, we can write v(r) as 




FIG. 4: Schematic picture of a region occupied by a loop, 
represented as a sphere with a radius of s y . Monomers in the 
region y > enter into this domain along the s-axis, whereas 
segments in y<0 along the —x direction. 



where V is the quantity having the dimension of the vol- 
ume (see (4.17)). 

b. From Loops to Dangling Chains A loop becomes 
a dangling chain with the transition rate u by detaching 
its one end from the junction with which its both ends 
has stuck. The disconnecting process of the loop are influ- 
enced by collisions with dangling chains and other loops. 

Since the transitions from the active chain to the dan- 
gling chain is also the disconnecting process of one end 
from a junction, the breakage rate 14(7) for the loop when 
7 = should be the sam e as t he detaching rate for the ac- 
tive chain appeared in ( 3.30| ), i.e., 14(7 = 0) = 0q. Let us 
introduce the quantity u(j)/f3o (77^0) which stands for 
the probability for the loop to dissociate during the time 
duration 1/0q. Note that l/0o stands for the lifetime of 
loops when 7 = 0. The ratio uiff)/ (3q must be larger than 
1, since the collisions of other chains with the backbone 
of the loop affects the possibility for its one of the ends 
to disconnect from the junction. We can assume that the 
increment A of w(7)//?o, defined by 



"(7) 
00 



1. 



(3.32) 



is given by the number of segments colliding with the 
loop during its lifetime. The collisions take place when 
the segments enter into the region occupied by the loop. 
We assume that the radius of the region is given by its 
radius of gyration. 

Under the above assumption, A is estimated as follows. 
Let the center of mass of a loop be the origin. Denoting 
by s y the radius of gyration of the loop of y-component 
(see ( p. 27 )), the number of segments vf^(y)dy belonging 
to dangling chains and loops, which enter into a region 
y^y + dy (y > 0) and —z(y)^z(y) occupied by the loop 
(see Figure 0) in unit time, is given by 

"&\v)dV = (" d + 2 z(y) dy ■ y x (y) 

= 2(u d + v l )Ny/( Sv y-y* jy dy, (3.33) 



v(r) = v 6(r)V, 



(3.31) 



where we used z(y) = {s y ) 2 — y 2 and v x (y) — jy. 
v l (i — d, I) is the number of i-chains per unit volume. 
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Integrating (3.33) from to s y and multiplying by 2 to 
take account of the region y < 0, we have the number of 
segments entering into the domain occupied by the loop 
in unit time: 

4 (v d + ^)a 3 7V 5 / 2 7, (3.34) 



3 • 36 3 / 2 



where we substitut ed Sy given by (3.27) for the second 
equality. Dividing ( 3.34 ) by u(j), we obtain 



• I = t4t(^ W)a 3 7V 5/2 — 4-. (3.35) 



162 



Putting (3.35) into ( |3.32 ), and solving it for u(j), we 
finally get 



u(j)/(3 = — V ~ , , (3.36) 



where $ 8 — Nv l a? is the volume fraction of i-chains. In 
the experiments we are going to analyze, <& a is much 
smaller than <I> d and as will be shown later. In this 



case, $ d + Q e in (3.36) can be replaced by the volume 
fraction $ of all chains. Note that there is no r depen- 
dence in u(j), since the end-to-end length of the loop is 
zero. 



IV. OBSERVABLE QUANTITIES 



/3 <T(r,i)+W> d (r,t), (4.4) 

d M) ] 



K{t)-— r — — V 



M a {v, t) - (p + v(r))</, d (r, t) + u4>\v, t), (4.5) 



^(r,t)=«(r)^(r,t)-«^'(r,t). 



(4.6) 



We w ill g ive their derivations in appendix |b| The equa- 
tion |4f3) corresponds to the stoc hastic Liouville equa- 
tion^ whereas the equation ( |4.5| ) is equivalent to the 
Fokker-Planck equationrB 

In the long time limit t — ► oo, <f) l (r,t) becomes in- 
dependent of time (we denote it as 4 (r)). Letting 
Jj<^ l (r, t) =0 (i = a, d, I) in the above equations, we have 



V- Ur^r) =~M a (r)+pcf ) d (r), 



(4.7) 



. K N \ k B T 
K-— r-— V 



<j> d (v) 



= /? o a (r) - (p + «(r))0 d (r) + u 0'(r), (4.8) 
v(r)</> d (r) = u<//(r). (4.9) 

With the help of (Q, we can eliminate 4> l (r) in Q)S 



K N \ k B T_ 
K : — I r V 



«^(r) 



Cn I (n 

= M a {v)-pcj ) d {v). (4.10) 



A. Zero-Shear Viscosity 

Here, we derive the zero-shear viscosity 770. Since we 
treat the case in which the shear rate 7 is much smaller 
than fa, [3 {r) c an be regarded as @q. We use the basic 
equations (2.13) wr itte n i n ter ms of the E ulcria n des crip- 
tion. Substituting (3.3), ( 3.13 ) and ( 3.17 ) into ( 2.13; ), we 
have a set of equations for (f>j(r, t), </>^(r, t) and <p l f (r, t): 



d_ 

di 



<^(r,t)+V- 



k{t)v 



f> a fM 



-M a f {v,t)+ P (j> d f {T,t), 



^(M) + V- 



(4.1) 

/W/M) - (p + »(r))^(r,t)+^(r,t),(4.2) 



St 



<f> l f (r,t)=v(r)<f> d (r,t)-ucl) l f (r,t). 



(4.3) 



Taking average with respect to the stochastic process 
R(t) in above equations, we obtain 



d_ 

di 



0°(r,t)+V- 



re(t)r 



<TM) 



_Z. Number of Chains in the Absence of Flow 

The number of i-chains per unit volume in real space 
is given by 



v i = J dr <j) l {r) 



(4.11) 



According to (2.1), we obtain the equation for the 
number-conservation 



n = v a + v d + v\ 



(4.12) 



where n is the total number of chains in unit volume: 

n = [ dr 0(r). (4.13) 



In the current situation of the slow shear flow, v % is ex- 
panded in powers of 7: v % = i / o+^i|7|+- ■ ■■ Since the zeroth 
order term v l contributes to t)q as seen later, we now de- 
rive the expression for Uq, i.e., the number of i-chai ns in 
uni t vol ume at the equilibrium state. Integrating (4.7) 
or (4.10) for k = , we have 



= -/3 <+p4 



(4.14) 
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Integrating (4.9) for k=0, we obtain 
v V<p d (v = 0) = Povl 



(4.15) 



where <^ (r) means <f> l (r) at the equilibrium state, and we 
put u — Po since 7 = 0. Let us now introduce ^(r) by 



$(r) = ^V(r). 



(4.16) 



Note that because of (4.11), ip(r) is normalized as 
fdmp(r) = l. Defining V by 



V 



ip(r = 0) : 



(4.15) reduces to 



vqVq = (3 is l . 



(4.17) 



(4.18) 



Solving ( |4.14| ) and ( J4 . 1 8j ) with the help of (|4.12|) , we ob- 
tain 



(4.19) 
(4.20) 
(4.21) 



1 + /3 + ft, 



l + ^(l + -f) 



where we put x = r±, y = T2 and z = r§. Solving the 
above equations, we get root-mean-square of the end-to- 
end vector for the active and the dangling chain in the 
form 



No 2 



2 / \2- 

P \ 1 2p\ 



(xy) a = 

(y 2 ) a = 



+2( 7 r)^j +— 
Na 2 



„2\d 



3 

No 2 
3 

A^a 2 



2_ 

A 



A y At \ A , 



+ yr 1 + 



A 



d JVa 2 



l + 2(yr) 2 < 
P 



3 ^ 1 + A 



(y 



2\d 



Na 2 



(z 2 r = (z 2 ) 



2\d 



Na 2 



, (others)=0 



(4.31) 
(4.32) 
(4.33) 
(4.34) 

(4.35) 

(4.36) 
(4.37) 



2. Root- Mean- Square of End-to-End Vector 

Hereafter, we denote the expectation value of the phys- 
ical quantity A(r) with respect to the i-chain as 



(A(r)Y = 



JdrA(r)^(r) _ 1 



±-JdrA{v)4> l (v). (4.22) 



fdrft (r) 

Multiplying the equation ( [4.7j ) and ( 4.10 ) by r a r/3, and 



integrating them with respect to r, we have 

-K al (r 1 rp) a v a - (rar-fYn^v - 

= -Po(r a rp) a u a + p(r a r p ) d v d , (4.23) 



K ai {r^rp) d v d - {r a r 1 np 1 ) d u d + — — (r a rp) 



\d,,d 



2k B T 



S a ^ d = fo{r a r p ) a v a - p{r a r fi ) d u d . (4.24) 



Substituting (3.25) for k, we obtain 

-2^{xy) a u a = -A(^ 2 )V a + p(x 2 ) d v d , 
-i(y 2 ) a v a = -fi {xy) a v a + p{xy) d v d , 
Q = -My 2 ) a v a +P(y 2 ) d v\ 

-^V + ^L(xy) V - 2-y(xy) V 

(,N (,N 



(4.25) 
(4.26) 
(4.27) 



= (3 (x 2 ) a v a -p(x 2 ) d v d , (4.28) 
^ (xy) d v d - j(y 2 } d v d = p (xy) a u a -p(xy) <V, (4.29) 
_2^ + 2^v (y2)V = /3o(y2)va _ p(y2)V) (430) 

(,N (,N 



where r is given by ( 3.1 5| ) 



3. Zero-Shear Viscosity 
The shear stress caused by i-chains is given by 



4y = »Hyf(r)) i - 



(4.38) 



We can use the linear form for the tension f(r) of the 
chain, since, when 7 <C (3 , the end-to-end length r of 
the chain stays sufficiently short during the deformation 
compared with its contour length. Putting /(r) = K^r 
into (4.38|), we have 



K N v l (xy)\ 



(4.39) 



With the help of ( 4.32 ), we obtain the shear stress caused 
by active chains in the form 



a a xy = k B Tv a 



O( 7 2 ),(4.40) 



where Pq is given by ( |4~19| ). Then, the zero-shear viscos- 
ity is obtained as 



r£ = lim = k B Tv* 
7^0 7 



1 

A 



T 1 



P_ 

A 



.(4.41) 
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Similarly, with the help of ( 4.35| ), the shear stress caused 
by dangling chains becomes 



< = k B Tvtir + 0(7 2 ), (4.42) 



where v d is given by ( 1.2C ), and the zero-shear viscosity 
is obtained as 



V $=lim?f=k B Tv d T 1 + f- . (4.43) 
The total zero-shear viscosity is given by 770 = t]q + i] d . 



4- Remarks on the Relaxation Time 

In the experiments we are going to analyze]!] the relax- 
ation time r (~ 1 s. see the next section for its esti- 
mation) given by ( 3.15 ) is always much shorter than the 
specific time 10 -3 ^ 10 _1 s) of the shear flow. In 



this case, as we can see from ( 4.4C ) and ( 4.42 ), the effect 
of r becomes very week for 7 smaller than l/r(~ 10 5 s). 
Consequently, we will adopt the limit t — > from now on. 
In this limit, the dangling chain has always stress-free 
conformation which is regarded as the Gaussian.Lj Then 
(f> d (r) is written, even for finite 7 (see ( 4.1 6| ) ) , as 



where 



.3/2 



^(r) 



2nNa 2 



exp 



3r 2 
2Na 2 



The total zero-shear viscosity becomes at r — > 



Vo = Vq- 



(4.44) 



(4.45) 



(4.46) 



B. Non-Newtonian Viscosity 

Here, we investigate the shear viscosity for the shear 
rate including larger than XJBq. To this end, we ana- 
lyze the basic equations-, (2T2) written in terms of the 
Lagrangian description,E3 since this description connects 
</> a (r,i) (and <j> l {r,t)) with ip(r). 



1. Equation for Active Chains 

The e quati on for the active chains is given by putting 
i = a in ([2.12]) as 



J a (t,t )r(r a (t;r ,t ),t) 



exp 



- / dt' /3(r a (i';r ,i )) 

Jt„ 



+p f dt' exp - f dt"l3{r a {t";r ,t Q )) 

xJ a {t',t Q )v d {t')^ a {t';r ,t Q )), (4.47) 
where we used the relation 

4> d (r,t) = v d {t)i>{r), (4.48) 



ip(r) is given by ( |4.45| ). Recall that r°(i;r ,io) in ( |4.47|) 
is given by ( |3.4| ) . In order to obtain the equation for the 
expectation value of A(r) wit h resp ect to the active chain, 
let us multiply the equation (4.47) by A(r a (t; ro, to)) and 
integrate it with respect to r : 



dv A{v)4> a {Y,t) 



dr exp 



dt'(3(r a (t'-r,t Q )) 



A(r a (t;r,t ))r(r,t ) 



pj dt' J drexp 



dt"f3(r a (t";r,t')) 



xA(r a (t;r,t'))v d {t')yj(r 



(4.49) 



In deriving (4.49), we changed the integrati on va riable 
from r a (t; Tq, to) to r in the left-hand side of ( 4.49| ), and 
from r a (t'; ro, to) to r in the second term of the right-hand 
side. Under the steady flow, due to X(t;t') = A(i — t'), 
(|4.49D reduces to 



dr A(r)(j) a (r,t) 



dr exp 



dt'f3(r a (t'-rj) 



A(r a (*-*o;r))</- Q (r,io) 



rt-t( 

-p I dt' I c?rexp 



dt"(3{r a {t"-r)) 



xA(r a (t';r))u d (t-t')if)(r), 
dr e- Mt - to) A{r a {t - t ;r))(t) a (r,to) 

D' r «(t-t ) 

+ p[ dt' f dr e-^'A^t^r))^^-*')^!-)- (4-50) 

JO JD',{t') 



where 



and 



r a {t-r) =A(t)r 



(4.51) 



D' r , (t) = {r I r a {t'- r) < r* (W < t)}. (4.52) 

A steady state can be obtained by taking the limit 
t — ► 00 and to — > -co. Clearly, the first term in the 
right-hand side of ( 4.50| ) becomes in these limit. As 
for the second term, because of £— »oo, v d (t — t') can be 
regarded as the value v d in the steady state for t' < 00. 
On the other hand, for t' —* 00, the integrand becomes 
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exponentially small due to the factor e~^ ot . Therefore, 
v d (t — t') is al ways regarded as v d in the limit of steady 
state. Hence, ( |4. 22| ) for i = a becomes 



d roo 



(A(r)) a = p— dt dr e- 0ot A(r a (t; r))V>(r). (4.53) 
va h JD',{t) 



Hereafter, we will take 



D r ,(t) = {r\ r a (t;r)<r*} 



(4.54) 



instead of D' r , (t) as the region of the integration in ( f .53] ) 
as far as the steady shear flow is concerned. Then the 
contribution to the integration (4.53) from the region out- 
side D' r »{t) (i.e., r — r a (t — 0; r) >r*) emerges. Recall that 
the Gau ssian distribution function ip{r) in the integrand 
of (4.53) mainly weighs the contribution from r smaller 
than y/Na. Since it turns out that r* is close to Na, 
the contribution from the region r > r*(~ Na 3> vNa) 
is almost negligible. Therefore, we can substitute safely 
the region D r *(t) for D' r ,(t). 



2. Equation for Loops 

The equation for the loops is written as 
j\t, t )<f> l (r l (t; t , r ), t) = e^*"*" V/(r , to) 
+ [ dt'e-^-^ J l {t\to)v d {t') 



y.v{r\t'; t , r ))^(r z (i'; t , r ), *')• (4.55) 

Following the same procedure as the case of active chains, 
we have the equation for the expectation value of A(r) 
with respect to the loop: 

J drA{v)(j) l {v) = J dtjdr e~ ut v{r)A(r l (t; r))^(r) 



V 1 u 



Vip(0)A(Q) 



V u 



(4.56) 



where we used V given by (4.17) 



3. Number of Chains 



Putting A(r) = 1 in (4.53), we have 

v a =pv d f dtj dre-' 3ot iP(r) = v d ^, (4.57) 

J0 JD r ,(t) C 



where we introduced 
1 



dt dr e"' 3ot V'(r). (4.58) 
C Jo Jd t , (t) 



Note that when 7 —> 0, we find £ — > [3q. Similarly, by 
putting A(r) = l in ( 4.56| ), we have 



4V0 



(4.59) 



By solving (4.58) and (4.5E) with the help of (4.12), the 



numbers of chains per volume for each type becomes 

n 



l + £(l + &y 



1 + H + M 



1 + ^(1 + !)' 

^. Shear Viscosity 
The shear stress caused by active chains is 



(4.60) 
(4.61) 
(4.62) 



pv 



a 



dt 







D r .(t) 



r a (t; r) 



where i/ d is given by (4.61) and r a (t, r) is given by (4.51). 
Recalling that the shear stress caused by dangling chains 
is zero, we finally get the formula of the total shear vis- 
cosity in the form 



7 

Tanaka-Edwards Limit 



(4.64) 



Tanaka and EdwardsB calculated the steady shear vis- 
cosity for the physical gel, composed of active and dan- 
gling chains (loops are absent) , with the assumptions that 
1) the breakage rate is given by f3(r)=j3 a +/3ir 2 (/3 and 
Pi are constant), and 2) even the active chains can be de- 
scribed by Gaussian chain model with the Hookean force 
f^{r) = KNr. The assumption 1) indicates that there are 
finite populations of active chains nearly equal to and/or 
even longer than the contour length I, i.e., there is no 
cutoff length r* (< I) in the TE model. The Gaussian 
chain model, proposed in the assumption 2), means that 
there are only shorter active chains (r<g^l) satisfying the 
Hookean force assumption. In this sense, the above two 
assumptions contradict each other. 

This conflict between the two assumptions 1) and 2) 
can be resolved by introducing the cutoff length r* much 
smaller than I. However, in general, there seems no rea- 
son to restrict the force to be linear. That is the reason 



why we adopt the force (3.1) given by the random-flight 



13 



model with the cutoff r* (0 < r* < Z). Hereafter, we call 
the case where r* <C Z and the TE limit within our 

model, since, in this case, the Hookean force is valid and 
loops are absent (u l =0). 



V. ANALYSES OF EXPERIMENTS 



No w, w e analyze, with the help of the formulae (4.64) 
and ( 4.63 ) derived theoretically in this paper, the stead 
shear viscosity observed experimentally by Jenkins et al. 
for HEUR having hexadecyl end groups in water. 



B. Estimate of the Relaxation Time 

We mentioned in the prece ding section that t he re lax- 
ation time r given by ( 3.15| ) (or n given by ( 3.2S| )) is 
much smaller than the specific time of the flow. This is 
quantitatively shown for HEU R by p utting the values of 
a and N listed in Table | into ( |3.29[ ) with T = 298K and 
% = 0.89 x 10~ 3 Pa-s (the viscosity of pure water at298K, 
latm). The obtained relaxation time for each molecular 
weight is shown in Table ||. Note that they have the 
order of 10~ 5 s. On the other hand, as seen from Figure 
|lC| , the interesting phenomena occur at a characteristic 
time in the range 1CP 3 ~ 10 _1 s. This is the quantitative 
reason why we took the limit r (or t\ ) — > in the previous 
section. 



A. Determination of N, a and $ 

The contour length Z of HEUR is obtained if we know 
the fully extended length Zo of ethylene oxide (EO), the 
main component of HEUR, and the number n rep of EO 
units in HEUR, i.e., Z = Zoxn rep . Taking the value 1.54A 
for a carbon-carbon bond length, 1.43 A for a carbon- 
oxygen bond length, and 70.5° for the angle formed by 
two bonds (tetrahedral bonds) in polymer backbone, we 
have Z = (1.54 + 2 x 1.43) xcos(70.5°/2) = 3.6A. Letting 
Mo (=44) stands for the molecular weight of the EO unit, 
we have n rep = M n /Mo, where M n is the molecular weight 
of HEUR. The values of I are listed in the. third column 
of Table | for each M n of the experiment. H 

Knowing the valueSrOf the root-mean-square y/ (r 2 )o of 
the end-to-end lengtbEJ in addition to I, we can obtain 
the number N of segments per chain and the length a of 
each segment by the relations N = l 2 /(r 2 ) and a = l/N. 
Here, (• • -)o indicates the average taken in equilibrium 
state without any flow. The values of N and a are listed 
respectively in the fifth and sixth column of Table Q. We 
see that the length a of the segment is almost indepen- 
dent of molecular weight M n (or N). Since the fully ex- 
tended length Zo of the EO unit is 3.6A, we can conclude 
that one segment contains 3^4 units. 

The volume fraction of polymers is defined by $ = 
Nna 3 with n being the number density of polymers. 
Since the number density is related to the mass density p 
of polymers by n — pNj^/M n , the volume fraction reduces 
to $ = 10 4 cNANa 3 /M n where Na is the Avogadro con- 
stant. Note that the mass density p is related to the poly- 
mer concentration c expressed by the weight percentage 
by (0 = 10 Xc (g/m 3 ). For each M n , the volume fraction $ 
of the c= lwt% HEUR aqueous solution is listed in Table 
H They are almost the same as they should be. We adopt 
the mean value $ = 8.0 x 10 -2 of them in the following. 
Note that the overlap threshold of polymers, estimated 
as $* ~ Na 3 /(r 2 f Q /2 ~ N' 1 ' 2 , is roughly 7x 10~ 2 when 
N = 227, for instance. It indicates that the volume frac- 
tion of polymers in this system, <!> = 8.0x 10 -2 , is on the 
order of the overlap threshold. 



C. Adjustable Parameters 



In our m odel, the shear viscosity rj is given by (4.64) 
with (4.63). We see that there are four parameters in 
(4.63), i.e., (3q, p, vq and r* (u is represented by (3q). 
Making use of the zero-shear viscosity obtained from the 
experiments!] we can eliminate one of three parameters 
Po, P, vq through the relation ( 4.46] ). For instance, p is 
represented by (3q and vq as 



P 



o + «o 



nk B T/r) f3 Q 



or, written in terms of (3q and vo/p as 

A) 



p 



nk B T/rj f3o - 1 - (vo/p) 



(5.1) 



(5.2) 



Hereafter, we adopt (3q and vo/p (and r*) as adjustable 
parameters. The ratio v /p is a measure representing 
which channels of the transitions are favorable, i.e., one is 
to the active, and the other to the loop from the dangling 
chain. 



D. Shear Viscosity 

In order to clarify the problems to be resolved, let us 
see the difference between the viscosity obtained from the 
TE limit (see Figure [l(] in appendix |(]) and the exper- 
imentally observed one by subtracting the former from 
the latter for each molecular weight, and by translating 
along the abscissa so that its position of the (first) peak 
coincides with others. As can be seen in Figure |[ there 
appear two peaks. The set of peaks at lower shear rates 
((i), in Figure ^|) reflects the shear-thickening behavior. 
The height of the peaks becomes low with increasing M n . 
It seems that each peak consists of two elements, i.e., one 
is high and narrow with strong M„-dependence, whereas 
the other is low and broad with weak M„-dependence. 
In the set of peaks at higher shear rates, i.e., at shear- 
thinning region ((h), in Figure o), the height of the peaks 
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M n 


Tlrep 

= M n /M 


i(A) 

— l()7l re p 


V<r 2 )o (A) 


N 

= W>o 


«(A) 
= l/N 


$ 

= W l cN A Na 3 /M n 


34,200 


777 


2,800 


186 


227 


12.3 


7.4xl0" 2 


51,000 


1,160 


4,180 


231 


327 


12.8 


8.1xl0~ 2 


67,600 


1,540 


5,540 


265 


437 


12.7 


8.0xl0" 2 


84,300 


1,920 


6,910 


301 


527 


13.1 


8.5xl0 -2 



TABLE I: Properties of HEUR in water. M n : number-average molecular weightpEj n rep : number of the EO units, I: contour 
length, -J (r 2 )o: root-mean-square of the end-to-end length] 2 !! N: number of segments, a: length of each segment and $: volume 
fraction for c = lwt% solution (mean value= 8.0 x 10" 2 ). Na = 6.02 x 10~ 23 is the Avogadro constant. For EO, Mo = 44 and 

; =3.6A. 



M n 


Tl(s) 


34,200 


1.3x10" 


5 


51,000 


3.1x10" 


5 


67,600 


5.4x10" 


-5 


84,300 


8.6x10" 


5 



TABLE II: Rouse relaxation time of HEUR obtained from 



(3.29) 
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FIG. 5: The difference between theoretically (TE limit) and 
experimentally obtained shear viscosity. The theoretical re- 
sults (lines in Figure [l^ in appendix Jo) are subtracted from 
experimental data (dots in Figure |l"o| ) for each molecular 
weight M n (N). To clarify the M n dependence of the height 
of peaks, each data is translated so that the positions of the 
first peaks (i) coincide with each other. 



has opposite M n dependence compared to the previous 
set of peaks. We will study mainly in this paper the first 
peak in the following. 

Let us here estimate r* from the viewpoint of bond- 
ing energy among hydrocarbon end groups. We considp; 
the potential barrier with height W around a junction.Ea 
An active chain dissociates from the junction when its 
one end climbs over the potential barrier W by the force 
f(r). Since the effective range of the hydrophobic inter- 
action, a, is sufficiently shorter than r* , we can put a = a. 



Then, the cut-off length r* is roughly estimated by the 
condition 



W = 



dr' -f(r') ~ f(r*)a. 



(5.3) 



The right-hand side of (|T^) represents the work done 
by the tension f(r) to detach the end from the junc- 
tion when r = r* . Assuming that r* < I, we have from 
Q W/k B T~(l~r* therefore, 7* /l~l-k B T/W. 
According to Annable et alp W = 28k B T (T = 298K) 
for HEUR of M w ~ 35, 000 having hexadecyl end groups 
(called C16/35 in the literature). In this case, r* is esti- 
mated as r* ~0.96L 



1 . Effects of Loops 

The number of each type of chains are plotted in Fig- 
ure | for r* = 0.97Z, A^Ss^ 1 and w /p = 245 (A^ = 227). 
With increasing the shear rate, the number v l of loops 
decreases and the number v d of dangling chains increases, 
since the transition probability u from a loop to a dan- 
gling chain increases with shear rate (see (3.36)). Due to 
the constant transition rate p from a dangling to an ac- 
tive chain, the number v a of active chains also increases 
with the same rate as the dangling chains at lower shear 
rates. That contributes to the gradual enhancement of 
the viscosity at lower shear rates as shown below. The 
enhancement of the number of active chains at the sJaear- 
thickening region has been observed by Tarn et al.U For 
higher shear rates, the population of active chains whose 
end-to-end length reaches r* becomes large, which causes 
the decrease in v a leading to the shear-thinning. The 
numbers of active and dangling chains and the number 
of loops in these figure s ind icate that the assumption 
v a <^.n d and v a <g;n' in (3.36) is reasonable (it also holds 
for other values of N). 

The steady shear viscosity for r* =0.97Z, /?o = 8s -1 and 
i; /p = 245 (iV = 227) is shown in Figure 0. We see that 
the calculated viscosity coincides well with experimen- 
tal data for the broad shear rate. Note that the value 
r* = 0.97Z is quite close to the value 0.96Z derived from 
the consideration based on the potential barrier between 
a junction and a hydrocarbon end group. In the figure, 
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FIG. 6: The number of active chains, dangling chains and 
loops plotted against the shear rate for r*/7 = 0.97, flo — 8s" 1 
and «o/p = 245 (TV = 227) corresponding to Figure m 
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FIG. 7: The steady shear viscosity plotted against the shear 
rate for TV = 227 (M„ = 34, 200). The solid line represents 
the theoretical result, whereas the dots represent the exper- 
imental results observed for 770 = 3.0Poise, T = 298K and 
c=lwt%£l 



we see the gradual enhancement of the viscosity at lower 
shear rates. This is because the number of active chains 
increases with shear rate in this region as shown in Fig- 
ure |6[ In addition, there is the effect of stretched active 
chains. It adds the sharp enhancement in the viscosity 
at the shear rate where the end-to-end length of the ac- 
tive chain reaches r*. When loops are absent, we had 
to adopt r* very close to I (see appendix |(]).E3 However, 
in the present case, since the increase in the number of 
active chains at lower shear rates promotes the viscosity 
rising, we can adopt proper value for r* . 

One might suppose that the observed shear viscosity 
can be fit without considering the high extension of active 
chains. However, it is impossible. When r* = 0.5/, for 
instance, we cannot heighten the peak, up to t he h eigh t 
of observed one with keeping the relation (5.1) or (5.2). 
Therefore, in the system we are dealing with, the above 



FIG. 8: The steady shear viscosity plotted against the shear 
rate. The lines represent theoretical results obtained for each 
M„ listed in Table |. The dots represent tine experimental 
results observed for T = 298K and c=lwt%£l The zero-shear 
viscosity r] is 3.0Poise for M n = 34, 200, 0.94Poise for M n = 
51,000,p0.53Poise for Mn = 67,600 and 0.26Poise for M n = 
84, 300-tl 



two mechanisms (i.e., the enhancement of the number of 
active chains and the high-extension of active chains) are 
indispensable. 

The steady shear viscosity are calculated as shown in 
Figure || with optimal values of r* , (3q and vo/p (given 
inside the box in the figure) for each molecular weight 
listed in Table D, We see from the figure that the theo- 
retical curve for each molecular weight M n fits very well 
with the first peak (although when M n = 84, 300, it does 
not fit well at higher shear rate). Note that r* does not 
depend on N, except for TV = 527. 

For each vq/p, the values of p and vq are given via 
J$), e.g., p = 0.14s" 1 , v = 35s" 1 for v /p = 245 (N = 
227), p = 0.1s" 1 , v Q = 17s- 1 for v /p = 165 (N = 327) 
and p= 0.065s" 1 , wo^e.Ss" 1 for v /p=105 (N = 437). 
Assuming that p and vq obey, respectively, the scaling low 
with respect to N, i.e., p cx iV" e and v cx N~**, we see 
from these data that e~ 1.2 and £~2.5 by the method 
of least squares. The dependence of p on N may be 
interpreted as follows. Since the polymer concentration 
is fixed, the number of hydrophobic end groups within 
the system decreases as the molecular weight increases. 
It diminishes the chance for a dangling chain to become 
an active chain. For example, when the molecular weight 
is doubled, the number of end groups within the system 
is reduced by half. This suggests p cx ./V -1 giving e = 1, 
hence the value of the exponent e~ 1.2, obtained by the 
present analysis, is fairly reasonable. The dependence 
of vq on N is also fairly reasonable since the probability 
for a self-avoiding chaia to form a loop is proportional 
to TV" 2 giving £ = 2a Recall, however, that we used 
just three data to determine the exponents. Obviously, 
we need more data in order to discuss further details of 
the N dependence of p and vq . It is desirable that some 
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experiments will be conducted to provide enough number 
of data for different molecular weights which may support 
a more precise analysis on the N dependence. 



VI. SUMMARY AND DISCUSSION 

We developed the network model for physical gel com- 
posed of unentangled linear chains with associative func- 
tional groups at both their ends. The presence of three 
types of chains was assumed in the network system: ac- 
tive chains with both their ends attached to the net- 
work responsible for the stress of the system, dangling 
chains with one of their two ends sticking to the network 
and loops with both their ends connecting to a junction. 
We investigated the transition rates between these chains 
how they depend on the end-to-end vector r and on the 
shear rate 7. The introduction of cutoff length r* was 
shown to be essential for the treatment of stretched ac- 
tive chains. It was also revealed that the transition from 
a loop to a dangling chain is promoted by the shear flow 
because of the collision with dangling chains or other 
loops. 

On the basis of the above network model, we showed 
that the origin of shear-thickening behavior is ascribed 
to 1) the dissociation of loops enhanced by shear flow, 
and 2) the high extension of active chains promoted by 
the flow. The portion of low and broad peak with weak 
TV-dependence originates from the mechanism 1), and 
the portion of high and narrow peak with strong N- 
dependence stems from the mechanism 2). 

It is revealed that the first peak in Figure || can be well 
explained by these two mechanisms. In Figure ^|, we put 
the difference between theoretically obtained shear vis- 
cosity (Figure |J) and the observed shear viscosity. The 
disappearance of the first peak in the difference declares 
that we succeeded to explain the shear-thickening phe- 
nomenon by the present theoretical analysis. 

In our analysis, we took the limit r — > for the re- 
laxation time of dangling chains, since the characteristic 
time of the shear flow is around 10 -3 ~ 10 _1 s, while r 
is estimated as 10 -5 ~ 10" 4 s (see Table [nj). Within the 
limit r — > 0, an active chain which detaches its pie end 
from the network relaxes to the Gaussian chainscj in an 
instant, and never rejoins to the network during the re- 
laxation. If the finite relaxation time are taken into ac- 
count, on the contrary, the detached chain has the possi- 
bility to recombine with thej-aeLwork before it totally re- 
laxes to the Gaussian chain.E3c!J Note that r is larger for 



larger N (see (3.15)), implying that the detached chain 
with larger N has higher possibility to recombine with 
the other junction during its relaxation. It indicates that 
the chains with larger molecular weight contribute more 
to the viscosity around the shear rate where characteris- 
tic time of the flow is the same extent as that of dangling 
chains. We suppose that this mechanism may explain the 
second peak (denoted by (ii) in Figure || and Figure |9|), 
since it appears at higher shear rates (7 > 10 3 s _1 ), and 
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FIG. 9: The difference between experimentally and theoreti- 
cally obtained shear viscosity. The theoretical results (Figure 
^) are subtracted from experimental data. To clarify the TV 
dependence of the height of peaks, each data is translated so 
that the positions of the observed peaks coincide with each 
other. See also Figure H for a reference. 



the larger the molecular weight is, the taller the height of 
the peaks becomes. To include these mechanisms into the 
present formula is one of the interesting future problems 
and will be reported elsewhere. 

Finally, let us comment on the idea of 'shear-induced 
transfer of intramolecular to intermolecular associations', 
proposed and investigated by Witten and CohenEJ in 
qualitative ways, and— developed by Ballard et al. in 
quantitative mannersJiil They considered solutions of as- 
sociative polymers having some functional groups along 
chains. When a flow is not added to the solution, each 
chain is in the random coil state, and hence, has a num- 
ber of associations among functional groups inside the 
chain (intramolecular associations). On the other hand, 
when sufficiently large velocity gradient is added to the 
solution, the chain extends due to viscous interactions 
with the solvent, thereby breaking there intramolecular 
associations. This causes associations among functional 
groups belonging to different chains (intermolecular as- 
sociations), and enhances the viscosity. If the number 
of functional groups is two, and these functional groups 
are attached at each end of a chain, the intramolecular 
association makes a loop. In this case, shear-thickening 
is caused by 'shear-induced transfer of loops to active 
chains via dangling chains' in their interpretation. 

One might say that this process is the same as the pro- 
cess proposed in the present paper. However, the origin 
to open a loop is quite different, i.e., as shown in section 
III , the dissociation of a loop is caused by collisions with 
other chains, in the present model. Indeed, viscous inter- 
actions with solvent affect the conformation of loops (see 
( 3.26 )), and enhance the possibility for a loop to be a dan- 
gling chain, however, this works effectively only at higher 
shear rates (7 > 1/Vi ~ 10 4 s _1 ). Therefore, it seems that 
the dissociation of loops caused by viscous interactions 
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do not contribute to the shear-thickening behavior at the 
moderate shear rates, although they may influence the 
second peak in Figure |[ To include the mechanismeS 
into our model is also interesting future problems. 



APPENDIX A: RADIUS OF GYRATION OF A 
LOOP 

We use the sine trans form for x(n, t) so as to satisfy 
the boundary condition (3.22), i.e., 



x(rc,t) = 2^x p (t)sin^— p 



(Al) 



We also use the sine transform for H(n,t), i.e. 



R(n,t) = 2^R p (i) sinj 

P =i 



N J 



Due to ( 3.20 ) and ( 3.21 ), the average and the variance of 
R p (i) are given by 



((R pa (t))) = 0, 



(A3) 



((R pa (t)R qP (s))) = ^Ls aP 5 m 5{t - s), (A4) 



respectively. Substituting ( |A1| ) and ( |A2| ) into (3.19), we 
obtain equations for x p (<) as follows: 



dt 



k(t)-—)x p (t) + ^-R p (t), (A5) 



where we put relaxation time of 'mode' p as 

_ CiN 2 1 _ CiA 2 a 2 1 
Tp ~ ~K^ 2 ~~p 2 ~ ~ 3TT 2 k B Tp*' 



(A6) 



Integrating (A5), from time — oo to t, we have 



Making use of (|A4|), the correlation functions of x pa {t) 
are written as 



({x P a(ti)x q p(t 2 ))) = 5 pg-^ 



kgT / -min ( i l>*2) 



dt / e _ (tl+ta _2t')/r, 



xJ2^ 1 (ti,t / )\ l 3 1 (t 2 ,t'). (A8) 

7=1 

For the stead y shear flow represented by the deformation 
tensor ( 3.24 ), we have 



Xpx(h)x qx (t 2 )}} = S, 



No 2 1 



\tl—t 2 \/T p 



^ + \r p ^\t-s\+ l -(T^) 2 Y (A9) 



((Xpy(h)Xgy(t 2 ))) = ((x pz (ti)x qz (t 2 ))) 

No 2 1 



= 5, 



((x px (h)x qy (t 2 ))} = 5. 



No 2 1 



*1 — *2 |/Tp 



|*x — t2\/r p 



(A10) 



X7((*i-t2)0(tx-*a) + fj, (AH) 
((x px (h)x gz (t 2 ))) = ((x py (h)x gz (t 2 ))) = 0. (A12) 
Putting 



x G (i) = — / dn x(n,f) = - V 

p:odd integer 



X p (t) 



(A13) 



(A2) into ( |3.23D , we get 



s a S(3 = 2^2((x pa (t)x p p(t))) 

p=l 



1G 



E 



(*)^/j(*))> 



>.g:odd integer 



pq 



(A14) 



Su bstitu ting (|A9|) - flAl2| ) into ( [ALU) , we obtain ( j3~26[ ) 

~ (FH). 



*p(t) = 7-[ dt' e-^'^XMKpit'). (A7) dn( t) 

S 1 J —oo 



APPENDIX B: DERIVATION OF THE 
FOKKER-PLANCK EQUATION 

Since <j>\ (r, t) is a stochastic variable, its differentiation 
with respect to time is n ot w ell-d efined. So, we should re- 
place the equations (1.1) ~ (4.3) by the forms of stochas- 
tic differential equations. As for (fh^), we have 

#/(r,t) = dfi(r,*)o<^(r,*) + # / (r,*)dt, (Bl) 

where Q(r,t)dt is the stochastic operator defined by 

1 



-V- 



k(t) - — - ) rdt + -^dw(t) 
(n I (n 



, (B2) 



and 



(r, t) = M a f (r, t) - (p+v{v))4> d f (r, t) +ucf> l f (r, t) . (B3) 



w(i) in (B2) is called the Wiener process. The correla- 
tions of its increment dw(t) are given by 

((dw a (t)))=0, (B4) 
((dw a (t)dw (s))) = 2( N k B T6 al3 S(t - s)dtds. (B5) 

Note that the Wiener process w(i) is related with the 
Gaussian white process R(f) through 

w(t) = f dt' K(t'). (B6) 
Jo 
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Using the fact 

dw a (t)dwp(t) — 2(^NkBTd a /3dt, 

we obtain 

1 

'Cat 



(B7) 



dti(r, i)dfi(r, i) = ^ a \7 p^dw a (t)dw (t) 



a, @ 

V^dt, 



(B8) 



where we neglected the terms of order higher t han dt for 
the first equality. Note that the product rule (B7) holds 
in the sense of stochastic convergences The relation (|B8| ) 
is the content of the fluctuation- dissipation theorem of the 
second kind. Hereafter, we will omit the argument r in 
(jA(jr,t), \&/(r, t) and fl(r, t) for simplicity. 

In (Bl), there appears the product of two stochastic 
variables, i.e., 4>At) and dw(t). There are two Types of 
such product: one is the product of the Ito type23 which 
is defined by 

<P d f (t) ■ dw(t) = <f>j(t)[w(t+dt)-w(t)], (B9) 

and the other is that of the Stratonovich type@ defined 
by 



4> d f {t) 



d> d (t+dt)+^ d (t) 

o dw(t) = — - [w(i- 



■cft)-w(f)].(B10) 



Since <fi d (t) is independent of the future stochastic pro- 
cesses w(t + dt) — w(i), the Ito product has a notable 
property 

• rfw(t))) = Uj(t))){{dw(t))) = 0, (Bll) 



where we used ( |B4| ) for the second equality. This prop- 
erty o ften sim plifies our calculations. From the defini- 
tions (B9) and ( |B10| ), we have the formula connecting the 
Stratonovich product with the Ito product in the form 



4>j(t) odw(t) = (f>j(t) • rfw(t)- 



~d4j(t)dw(t). (B12) 

As stated in subsection II D| , we regard the stochastic 
product appeared in (Bl) as that of the Stratonovich 
type 



Applying the formula (B12) to dfi(t) o <j£(t) in (Bl). 



and neglecting the terms of order higher than dt, we have 
the stochastic equation of the Ito type 



d(f> d f (t) = df2(t) ■ 4>f(t)+V f (t)dt, 



(B13) 



where 

dO(t) = dn(t) 
= dfl(t) 



-dn(t)dft(t) 
V 2 ^dt 




-V-dw(t). 
(B14) 
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FIG. 10: The steady shear viscosity plotted against the shear 
rate (TE limit). The lines represent theoretical (TE limit) 
results obtained for each M n listed in Table | with r*/7 = 0.5. 
The dots represent tfee experimental results observed for T= 
298K and c= lwt%£l The zero-shear viscosity 770 is 3.0Poise 
for M n = 34,200, 0.94Poise for M n = 51,000, 0.53Poise for 
M„ = 67,600 and 0.26Poise for M n = 84, 300.1 



Takin g the random average of the stochastic equation 
(B13), we obtain the Fokker-Planck equation 



d<f) d {t) 



-V- 



-f(t)dt, 




k B T. 



5 JV 



cj> d {t) dt 



(B15) 



In 



where 4> l (t) — ((4 > l f (t)) ) as defined in subsection II D 
order to obtain ( B15|) fro m (B13), we used the property 
of the Ito product^ Bll ) . 



As for the equations (4.1) and (4.3), since they do not 
contain the stochastic process R(t) explicitly, we can ob- 
tain (4.4) and (4.6) in a st raight forw ard way by taking 
the random average of (hi) and (4.3), respectively. 



APPENDIX C: SHEAR VISCOSITY (WITHOUT 
LOOPS) 

1. TE Limit 

Here, we compare the steady shear viscosity observed 
by Jenkins et alB with that given by the TE limit (Fig- 
ure [n]) for each molecular weight listed in Table |[ As 
mentioned in subsection |lV Cj , the TE limit is attai ned b y 
adopting small r* and put ting vq — in equation ( f4.63| ) . 
Since p is given through (5.1) or (5.2), the number of 
adjustable parameters is two, i.e., /3q and r*. 

We adopt r* — 0.51 here, although the value is not 
so small enough being good foi^the Gaussian model or 
the Hookean force assumptions With this r* , the value 
of (3q is determined by adjusting the theoretical curve 
to experimental data at the shear-thinning region (the 
abscissa is scaled by /3q for given N). For example, the 
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FIG. 11: The number of active chains plotted against the 
shear rate for N = 221. The solid line stands for the number 
of active chains within the TE limit which corresponds to the 
solid line in Figure [j^. The dashed line represents the number 
of active chains corresponding to the solid line in Figure [l^. 



FIG. 12: The steady shear viscosity plotted against the shear 
rate (without loops). The lines represent theoretical results 
(«o = 0) obtained for r*/l = 0.5, 0.9, 0.99 and 0.99994 with 
fixed N = 227 and /3o = 7s _1 . The dots represent the ex- 
perimental results obsewed for M n = 34, 200, rjo = 3.0Poise, 
T = 298K and c=lwt%£l 



adjusted value of (3 is 70s" 1 for iV = 227 (M„ = 34, 200) 
giving us the lifetime 1//3q — 0.01s of active chains. It is 
much smaller than the observed relaxation tiflfie of HEUR 
aqueous solutions having the order of 0.1s .Ero The values 
of Po for other N are given inside the box in Figure |l(| 
We see from these figures that the shear viscosity de- 
creases monotonically as the shear rate increases, i.e., 
the shear-thickening behavior does not appear within the 
TE limit. With increasing the shear rate, the number 
of longer active chains increases. The chains which ex- 
ceed r* dissociate from the junction and become dangling 
chains. This causes the decrease in the number of active 
chains (see the solid line in Figure |ll|.) leading to the 
thinning in the shear viscosity. Since r* is not close to I 
within the TE limit, the tension of the active chain with 
end-to-end length r* is not so strong. Therefore, they 
can not generate the enhancement of the viscosity. That 
is the reason why the viscosity rise does not occur in the 
TE limit. 



2. Effects of High-Extension of Active Chains 
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FIG. 13: The steady shear viscosity plotted against the shear 
rate (without loops.) The lines represent theoretical results 
(vo = 0) obtained for each M„ listed in Table Q. The dots 
represent rthe experimental results observed for T = 298K and 
c= lwt%Jj The zero-shear viscosity r/o is 3.0 for M n = 34, 200, 
0.94 for MnfF 51,000, 0.53 for M„ = 67,600 and 0.26 for 
M n = 84,300.tl 



Let us now see the influence of the stretched active 
chains on the viscosity. To this end, we consider the 
case where loops are absent (vq = 0) in the wake of the 
preceding subsection. 

The r* dependence of the steady shear viscosity is 
shown in Figure [l| for N = 227 (M„ = 34, 200) and 
(3q = 7s~ 1 . From Figure 12, we see that a peak starts 



to appear when r* > 0.5Z, and the height of the peak be- 
comes higher with increasing r* . When r* — 0.99994Z, 
it coincides with the experimentally observed one. Note 
that the appeared peak is narrower than observed one. 
The value /3q = 7s _1 is determined so that the position 
of the peak coincides with observed one. The obtained 



l//3o = 0.14s 1 has the same order_^s observed relaxation 
time of HEUR aqueous solutions .Bel 

Figure [l3] shows the theoretically obtained steady 
shear viscosity with optimal values of r* and (3q (given 
inside the box in the figure) for each molecular weight 
listed in Table |, compared with experiments. We see 
from these figures, in common, that 1) the cutoff length 
r* is too close to the contour length I, and 2) the peak is 
too sharp. 

The number of active chains for N — 227 is shown in 
Figure |ll| (the dashed line) . As one can see, it monotoni- 
cally decreases with increasing the shear rate, the reason 
of which is the same as in the case of r* =0.51 (i.e., TE 
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limit). In order to enhance the viscosity against the de- 
creasing number of active chains, we are forced to choose 
r* very close to I, since r* < I causes the strong force 
between ends of active chains. That is why r* so near 



to I is required to fit the experimental data. The above 
problems 1) and 2) is resolved by taking account of the 
open process of loops. 
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the long time limit. \~\n 
This limit was adopted in s ome literatu res EJ lij 
Note that we can_nse ( |2.12j ) instead of ( |2.14| ) when r is 0. 
In the literature,E3 four distinct values of -J (r 2 )o at 30° C 
are reported for each molecular weight, which were ob- 
tained using four different theoretical relations connecting 
(r 2 )o and the intrinsic viscosity [n] . The values of {r 2 )o 
listed in Table | are mean values of these data. Although 
the steady shear viscosity we are going to analyze is ob- 
served at 25° C, since [77], and therefore yj (r 2 )o, do not 
change their values so much within the temperature range 
25~30°C (see Figures 7,8 and 9 in the literaturecJ), we can 
use yj (r 2 )o at 30° C to analyze the steady shear viscosity 
at 25° C. 

At equilibrium state, the conformation of a dangling chain 
is regarded as the Gaussian. 

Their discussions are not based on the transient network 
theory. 

That is, the shear-induced transfer of loops to active chains 
via dangling chains, by means of frictional force due to 
solvent. 

This reason is as follows. For smaller r* , the viscosity starts 
to decrease at lower shear rates, since the lifetime of an 
active chain is shorter effectively in this case. It means 
that, for smaller r* , we have to adopt larger /3o (smaller 
l//3o) so that the theoretical results fit to experimental 
data at the shear-thinning region. Therefore, in order to 
get larger l//3o close to the observed relaxation time, we 
have to choose larger r*. Although the value r* —0.51 gen- 
erates stronger force compared with Hookean force, it does 
not cause viscosity rising. This value seems to be the crit- 
ical one for generating the shear-thickening behavior (see 
Figure |l§). 



